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§ 1 - Introduction. 

The study of semi-classical orthogonal polynomials in the scalar case (i.e., the 
study of orthogonal polynomial whose associated functional satisfies a distribu- 
tional equation _D(u$) = u*&, where $, ^ are polynomials with deg \& > 1) was 
started by Shohat ([15]) in order to generalize the properties of classical orthogonal 
polynomials. 

Among others, in ([2, 9]), an approach to such polynomials taking into account 
the quasi-orthogonality of the derivatives of the polynomials sequence is given. We 
also mention the results of Maroni ([11, 12]) where an algebraic theory of semi- 
classical orthogonal polynomials is presented. The purpose of this theory is the 
characterization of semi-classical orthogonal polynomials by means of the quasi- 
orthogonality of their derivatives, the structure relation and the differo-differential 
equation of second order. 

In the last years, the study of matrix orthogonal polynomials attracted a great 
interest of the researchers, (see [4, 6, 10]). As for their differential properties it is 
known the result of Duran ([5]) who characterizes those matrix orthogonal polyno- 
mials (O P) satisfying a symmetric second order differential equation with polyno- 
mial coefficients. He proves that they are diagonal (up to a factor) with classical 
scalar O P in the diagonal. In spite of the variety of applications of matrix poly- 
nomials ([4, 5, 6]), there are not too many known families of semi-classical matrix 
O P out of the diagonal case. 

One way to study many families of matrix orthogonal polynomials is to extend 
the analysis started by Duran of differential properties of matrix orthogonal poly- 
nomials. A natural way to do this is to generalize the theory of semi-classical scalar 
O P to the matricial case. In order to do that, we start with a Pearson type equa- 
tion for matrix quasi-definite functionals. We obtain their characterization in terms 
of a (in general, non-symmetric) differo-differential equation for the related matrix 
orthogonal polynomials. In the way of the proof, we find another equivalences that 
generalize the scalar case. 

This paper has been organized as follows. In Section 3 we introduce and study 
the concept of quasi-orthogonality for matrix polynomials. 

The mean result of Section 4 is the caracterization of semi-classical functionals 
in terms of the quasi-orthogonality of the corresponding matrix polynomials and 
their derivatives. 

In Section 5 we prove for matrix O P with the above quasi-orthogonality proper- 
ties, the structure relation and the differo-differential equation, showing that they 
are equivalent to a Pearson-type equation. 

Finally, in Section 6, we illustrate the preceding with some examples and we find 
a way to construct non-diagonalizable semi-classical matrix functionals. 

§ 2 - Basic tools. 



We shall denote by P( m > the C( m < m Ucft-modulc 
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and by means of P<" 1 )' the C( m ' m )-right-module Hom(p( m \ C^" 1 '™)) . 

2.1 Definition. 

(i) The duality bracket is defined by 

(•,•): P^ m ' > x p(" 1 )' — ► c( m - m ) 

(P,u)-(P,u) :=u(P) 

(ii) Fork G N and u <= P (m) ' the linear functional ux k I G P (m) ' is defined by 

(P,ux k l) := (x fc P,u), 

wfiere J denotes the m x m identity matrix. 

A linear extension gives the right-product u ■ A G p(" 1 ) for u G p( m ) j 

A e P("), with A(x) = 

akX k , in the following way: 

k=0 

n 

(P,uA) = Y,( xkp > u ) ak - 

k=0 

(Hi) The inner product 

/ A . p(m) x p(m) _^ £.(m,m) 

(P,Q)^(P,Q) u := (P,uQ*) 

is defined for every m G P^" 1 ' , wfiere Q* denotes the trasposed conjugated of Q. 
Remark. 

The duality bracket verifies the following linear properties: 

(«iPi+a 2 P 2 , P\U\+(3 2 u 2 ) = 

= ai(Pi,ui)/?i + a 1 (P 1 ,u 2 )p2 + a 2 (P 2 ,u 1 )(3 1 + a 2 (P 2 , u 2 )(3 2 , 

for all ai,a 2 ,/3i,/? 2 G C( m ' m \ Pi,P 2 G P (m) and ui,u 2 G F (m) ', while the inner 
product is sesquilinear: 

(i) (a 1 P 1 + a 2 P 2 , fcQi + p 2 Q 2 ) u = a 1 {P 1 ,Q 1 )jl + +a 1 (P 1 , Q 2 )J^+ 
+a 2 (P 2 , Q^Jt + a 2 (P 2 , Q 2 )J* 2 , 
for all ai , a 2 , 0i , /3 2 G C( m < m ) , Pi , P 2 , Qi , Q 2 e P (m) ; 

ai) <p,q>; = <q,p>„ 

for all P,Q G P (m) . 

2.2 Definition. We denote by Ck ■= (x k I, u) the k-th moment with respect 
to u G P(" 1 ' . Given u G P^™* 1 with moments (G\) keN , we say that u is quasi- 
definite (non-singular) if det[(Cfc+j)™ Q ] 7^ 0, Vn > 0, were (Ck+jYl j=o is ^ e 
Hankcl-block matrix 

/Co Ci • • • C„ \ 
Ci C 2 • • • C n +i 

\ C n C n +i ■ ■ ■ C 2n J 



4 DIFFERENTIAL PROPERTIES OF MATRIX ORTHOGONAL POLYNOMIALS 



Remark. Given the sequence 

(A)fclo c & m - m) , there exists a unique 
u G such that (x k I, u) = Ck- This establishes an isomorphism between P^" 1 ) 

oo 

and the formal series with coefficients in C^ m ' m \ C ^ 

k=0 

2.3 Definition. Let u G P< m )'. We say that u is hcrmitian if C* k = C k , for all 
k > 0. 

2.4 Theorem. Let u G P^™* 1 be quasi-definite. Then, there exists a unique (up 
to left non-singular matrix factors ) sequence of left orthogonal polynomials (Pn) n>0 
with respect to u, that is 

(1) P„eP( m ), dgP n = n. 

(2) The leading coefficient of P n is non-singular. 

(3) (P n , P m ) u = K n S nm , where K n is non-singular. 
This sequence verifies a recurrence relation 

xP n (x) = a n P n+1 (x) + (i n P n {x) + 7„P„_i(a;), n > 0, 

p-i(x) = e, 

with a n , p n , 7„ G C^ m ' m \ a n , 7„ non-singular and we denote by 9, the zero matrix. 
Proof. Sec [4, 14]. oo 

Remarks. 

(i) If we consider the structure of right-modulus on P^" 1 ) (so, left-modulus on P( m ) ), 
we define in a similar way a sequence of right orthogonal polynomials. 

(ii) The inner product (P n ,x n l) u is non-singular for every quasi-definite u G P^™* 1 . 



§ 3 - Quasi-orthogonality. 

Like in the scalar case, the semi-classical character will be closely related to the 
idea of quasi-orthogonality. But the scalar quasi-orthogonality, when generalized 
to the matricial case, splits in two different concepts. Both of them will play an 
important role in the characterization of semi-classical matrix functionals. 

3.1 Definition. Let (Qn) n>0 be a sequence of matrix polynomials such that 
the leading coefficient of Q n is non-singular and dg Q n — n. 

1. - We will say that (Qn) n>0 is a sequence of left quasi-orthogonal matrix poly- 
nomials of order r with respect to v G p( m ) \ {0} if 

(a) (x k Q n ,v) = 6, 0<k<n-r-l, n>r + l. 
(bl) There exists n > r such that (x na ~ r Q no , v) ^ 9. 

2. - We will say that (Qn) n>0 is a sequence of regularly left quasi-orthogonal 

matrix polynomials of order r with respect to v G p( m ) \ {0} if 
(a) (x k Q n ,v) = 6, 0<k<n-r-l, n>r + l. 
(b2) There exists n > r such that (x na ~ r Q na ,v) is non-singular. 

Notice that if r = the previous definition is the classical orthogonality. Fur- 
thermore, for m = 0, both definitions yield the scalar quasi-orthogonality. 
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Remember that given w £ P^" 1 ' we will say that w — if {P, w) = 9, 
VP € P (m) . 

As we will see in the next proposition, when a sequence of quasi-orthogonal 
matrix polynomials is already orthogonal with respect to another functional, the 
conditions (bl), (b2) become stronger. 

3.2 Proposition. Let u e p( m ) be quasi-definite and let (Pn) n>0 be the co- 
rresponding sequence of left orthogonal matrix polynomials. Given v € P( m ) , the 
sequence (Pn) n>0 is quasi-orthogonal of order r with respect to v if and only if 

(a) (x k P n ,v) =6, 0<k<n-r-l, n>r + l. 
(b'l) (x n - r P n ,v) ^6, Vn > r. 

The sequence (Pn) n>0 is regularly quasi-orthogonal of order r with respect to v 
if and only if 

(a) (x k P n ,v) = 6, 0<k<n-r-l, n>r + l. 
(b'2) (P n x n ~ r , v) is non singular Vn > r. 

Proof. Let n > r + 1 be verifying (bl) or (b2). Taking into account the 
recurrence relation (1.1), we obtain 

7n+i (x"- r P n ,v) = -p n+ i (x n - r P n+1 ,v) -a n+ i (x n - r P n+2 ,v) + 

+ (x n - r+1 P n+1 ,v) = (x n+1 - r P n+1 ,v) 
and (bl) or (b2) holds for n + 1. In a similar way, 

ln (x n - 1 - r P n - 1 ,v) = (x n - r P n ,v) 

and (bl) or (b2) is satisfied for n — 1. oo 

In what follows, the following result will be useful. 

3.3 Lemma. Let u e p( m )' be quasi-definite and let (P n ) n>0 be the corres- 
ponding sequence of left orthogonal matrix polynomials. Then, given w <G p( m ) , 
w — if and only if there exists no such that 

(x k P n ,w)=6, 0<k<n, n>n . 

Proof. As a consequence of the recurrence relation (1.1), 

< x k P no -x,w >= 

= ln^[< x k+1 P no ,w > -a no < x k P no+1 ,w > -0 no < x k P na ,w >}=6 

for < k < no — 1. So, the hypothesis is true for n > no — 1, < k < n too, and 
as a consequence, for all n > 0, k > 0. That means < x k , w >— 9 for all k > and 
then, w = 0. oo 

Given a sequence of orthogonal matrix polynomials with respect to a functional 
u, its quasi-orthogonality with respect to another functional v can be characterized 
by a simple relation beetwen u and v. 

3.4 Theorem. Let u, v G pW SU ch that u is quasi-definite and let (P n ) n>0 be 
the corresponding sequence of left orthogonal matrix polynomials. Then, 
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(i) (Pn) n>0 is quasi-orthogonal of order r with respect to v if and only if there 
exists A 6 p( m ), with dgA = r, and such that v = uA. 

(ii) {P n ) n >o i s regularly quasi-orthogonal of order r with respect to v if and 
only if there exists A e p( m ) ; with dgA = r, and non-singular leading coefficient 
such that v = uA. 

Moreover, in any case, the matrix polynomial A is unique. 
Proof. 

We will do the proof just for (i) because for (ii) the arguments are analogous. 

r 

<= / We consider A(x) = ^ ajx\ with a r ^ 9. Then, for v = uA 

r 

{x k P n ,v)=Y,{x k+i Pn,u)a i . 

3=0 

For < fc < n — r — 1, we have 

(x k+ ip n ,u)=6, (0<j<r) 

and, for k = n — r, 

(x k P n ,v) = {x n P n ,u) a r ^6. 
So, (P n ) n >o 1S regularly quasi-orthogonal of order r with respect to v. 

r 

=> / We will prove that = v — ^^x^uaj has an unique solution in the 

3=0 

unknowns (aj)'j =0 . 

For all n > r fixed, the system of r + 1 equations 

' (x n - r P n ,v) = (x n P n ,u)a< r n \ 

< (2-1) 

k (x n P n , v) = (x n P n ,u) 4 n) + ■ • • + (x n+r P n , u) 4"), 

has a unique solution in the unknows (a^)J =0 . Notice that the first equation leads 
to a { r n) + 0. 

We will prove that (a^ )j =0 are independent of n. Keeping in mind the recur- 
rence relation (1.1), we have that the relation 

7„+i (x k P n ,v) +/3„+i (x k P n+1 ,v) +a n+1 (x k P n+2 ,v) = (x k+1 P n+1 ,v) (2.2) 

is verified Vfc > and V£> € p( m ) . The first equation of (2.1) leads to 

f (x n - r P n ,v) = (x n P n ,u)a^ 



("+!) 



and, from (2.2) 



(x n+1 - r P n+1 ,v) = (x n+1 P n+1 ,u)a 



7n+ i (x n - r P n ,v) - (x n+1 - r P n+1 ,v) = 9, 
7„+i (x n P n ,u) - (x n+1 P n+1 ,u) = 9. 



(2.3) 



(2.4) 
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Equations (2.3) and (2.4) imply = ai n+1 \ 

We supose o[ n ^ = a|™ +1 ' (I = m+1, . . . ,r, m < r, Vn > r), and we will prove 
that a? = dm +1 '. Taking into account (2.1) we have 



{x n - m P ni v) = {x n - m+i P n ,u)af ] 

j=m 

r 

(x n - m+1 P n+1 ,v) = (x n+1 - m+ iP n+1 ,u)a 



(n+l) 
3 



and from (2.2) 



' 7 n+i (^"- ro P„,v) +/3„+i (a: n - m P n+ i,i;> + a n+1 (x n - m P n+2 , v) = 
= (x n - m+1 P n+ll v), 



(2.5) 



(2.6) 



' 7«+i (a;"- ro ^P n ,n) + f3 n+1 {x n ~ m+j P n+1 ,u) + a n+ i (a;"- ro ^P n+ 2, u) = 

= (x n+1 - m+ ip n+u u). 



By substitution of (2.5) into (2.6) we obtain 



(2.7) 



7 „ +1 ^<x"-^P„, u )4" ) +/3„ +1 ^ (^-^P„ +1 , u )4" +1) + 

j=m j =m +l 



+ an+1 Y J (x n - m+3 P n+ 2,u)af + ^ = £ {x n+1 - m+j P n+ uu)af +2 \ 

3=m j=m+2 

Now, keeping in mind the relation (2.7) and by application of the induction 
hypothesis, we have 

7n+ i (x n P n ,u)ag> -(x n+1 P n+1 ,u) a^ +l) =e. 
Thus, by (2.4), we conclude = flt +1 '- 

Let us denote by (aj)j =0 the solutions of (2.1) for n > r, that we have proved 

r 

are independent n. If A(x) = ajx-* and io = u — uA, we have 



(x k P n , w) = (x k P n , v) - (x k+j P n , it) aj, < fc < n, n > r. 

When 0<A;<n — r — 1, taking into account the orthogonality with respect to 
the functional u as well as the quasi-ortogonality with respect to the functional v, 
we get (x k P n ,w) = 9. If n — r < k < n, then (x k P n ,w) — 9, because (aj)J_ 
are solutions of the equations (2.1). So, w = 0, according to Proposition 2.2. 
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To prove the uniqueness of the matrix polynomial A, let us suppose A, B G p( m ) 
such that v = uA = uB. Then, from (i), it must be r = dgA — c\gB since (P n )n>o 
can not be quasi-orthogonal of different order with respect to v. Hence, A(x) = 

r r 

a,jx\ B(x) = bjX^ with a r , b r ^ 9. From 
j=o j=o 

p i 

(x k P n , V)=J2 (x k+j Pn, U) a, = ]T {* k+1 Pn, u) bj, 
3=0 j=0 

and taking k = n — r, . . . , n, we have 

' {x n P n ,u)(a r -b r )=9, 
< 

k (x n+r P n , u) (a r -b r ) + --- + (x n P n , u) (oo - b ) = 9. 
Therefore, aj = bj, (j = 0, ■ • ■ , r). oo 



§ 4 - Semi-classical functionals. 

We consider the derivative operator on the space p( m ) as the linear operator 
D : P<" 1 )' -» P(" 1 )' such that (P, Du) = - (P',u) . From this and the definition of 
the right-product it is straightforward to prove that D(uA) = (Du)A + uA', for 
all u e P (m) ' and A G P (m) . 

4.1 Definition. Let u G P^ m ) be quasi-definite. We will say that u is semi- 
classical if there exist A,B G p( m \ with det A =/= 0, sucii that it is verified the 
distributional equation D (uA) = uB. We will also say that the corresponding 
sequence of left orthogonal matrix polynomials (P n ) n>0 is semi-classical. 

Remark. Given u G p(" 1 )', the linear functionals u t j : P (m) — > C, = 
1 , • • • , m) defined by Uij(P) = u(P)ij, VP G P( m ) are called the components of 

m 

u. Then, (uA)ij = Ui^Akj with the obvious definition for the multiplication 

fc=i 

of Uij by a scalar polynomial. Therefore, detA ^ is the minimal requirement 
to ensure that the previous definition will involve to all the components of the 
functional u. 

4.2 Lemma. Let u G P (m) ' such that D (uA) = uB. Then, for every C G P (m) , 

D (uAC) = u (AC + BC) 

Proof. It is just a consequence of the rule for the derivation of the right-product. 

oo 

The previous result implies that, for every u G P^ m ' , the set 
An = {A e P (m) /D{uA) =uB,Be P (m) } 
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is a right-ideal of P< m ). 

When dealing with scalar semi-classical functionals we arrive in this way to an 
ideal of P, which is necessarily generated by a unique (up to non-trivial factors) 
polynomial. This generator is used to classify the scalar semi-classical functionals, 
([11, 12]). 

Since in the matricial case, a right-ideal of p( m ) is not necessarily principal, 
we cannot use A u for the classification of semi-classical functionals. However, if 
u E P^ 7 ™) is semi-classical then there exists a E P \ {0} (where P = P^) such that 
al E A u - To see this just notice that if A E A u with detA ^ 0, then (detA)I E A u 
since (detA)/ = AA+ with A+ =adjA e P (m) . Therefore, for every semi-classical 
matrix functional u E P( m ) , the set 

K = {a E P/D(ual) = uB,Be P (m) } 

is a non-trivial ideal of P. We can use the "essentially" unique generator of this 
ideal to clasify the semi-classical matrix functionals similarly to the scalar case. 

4.3 Definition. Let u E p( m ) be semi-classical and let a E P \ {0} be a 
polynomial with smallest degree such that D(ual) = uB, B E p(" 1 ). Then, we say 
that u is of class s —max{p — 2, q — 1}, where p — dga and q = dgB. 

Remarks. 

(i) The preceding discussion shows that this definition is well done: it always 
exists such a polynomial a and the definition of class does not depend on the choice 
of a. 

(ii) When u E p( m ) is semi-classical there exists A E A u with dctA ^= 0, but it is 
possible for the leading coefficient of A to be singular. However, there always exists 
A E A u with non-singular leading coefficient. To see this just take A =(detA)I. 

The following theorem gives the first characterization of semi-classical matrix 
functionals. 

4.4 Theorem. Let u E p( m ) be quasi-definite and (Pn) n >o^ e ^ e correspon- 
ding sequence of left orthogonal matrix polynomials. Then, the following statements 
are equivalent: 

(i) u is semi-classical. 

(ii) There exists v E p( m ) \ {0} such that (Pn+i) n>0 is quasi-orthogonal with 
respect to v, and (P n ) n >o is regularly quasi-orthogonal with respect to v. 

Proof. 

(i) (ii) Notice that if u semi-classical then D (ual) — uB, a E P \ {0} 

according to previous comments. Therefore, 

((x k P n )',ual) + (x k P n ,uB) =0. 
p q 

Let a(x) = ^ ajx\ B(x) = ^ bjx^ with a p ^ and b q ^ 9 whenever B ^ 8. 

3=0 j=o 
With this notation the above equation becomes 



p i 
53 (kx k - 1+ 3P n , u) a, + J2 {* k+3 Pn, u) bj = - {x k P; i7 ual) . (3.1) 

3=0 3=0 
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The left hand side is equal to zero matrix if < fc < n — p and < k < 
n- (q+ 1). Let r = max{p - l,q}. Then, {x k P^,ual) = 0if0<fc<n — r — 1. 
For k — n — r, (3.1) becomes 

(n-r) (x n P n ,u) a r+1 + (x n P n ,u) b r = - (x n - r P' n ,ual) , 

and thus 

(x n - r P^ual) = - {x n P n ,u} [(n - r)a r+1 + b r ] , 

which is equal to zero matrix for at most one n. So, the sequence (Q„) n>0 with 
Q n = P r ' +1 is quasi-orthogonal with respect to ual, of order r — 1. 

Obviously, since al has non-singular leading coefficient, (P n ) n>0 is regularly 
quasi-orthogonal with respect to ual, of order p. 

(ii) => (i) If {P n ) n >o i s regularly quasi-orthogonal of order p with respect to 
v, there exists A G P^" 1 ^ with dgA = p and non-singular leading coefficient of A, 
such that v = uA. Notice that if the leading coefficient of A is non-singular, then 
Act A ^ 0. 

Let w = D(uA). It is, (P, w) = - (P', uA) and for P = x k P n we have 

{x k P ni w) = -k{x k - l P ni uA) - {x k P' nl uA). 

From the quasi-orthogonality, the right hand side vanishes if < fc < n — p and 
0<fc<n — s — 2, where s is the order of quasi-orthogonality of (-F^+i) >0 - So, 
{P n ) n >o is quasi-orthogonal with respect to the functional w with order at most 
max{p — 1, s + 1}. Thus, there exists B E p( m ) such that dg-B < max {p — 1, s + 1} 
and w — uB. oo. 

Remark. Notice that if (-P,'+i)„ >0 is quasi-orthogonal of order s with respect 
to ual, then D{uat) = uB where dga = p, dgB = q and this implies that s = 
max {p — 2, q — 1}. 

§ 5 - Structure relation and differo-difFerential equation. 
5.1 Theorem. (Structure relation) 

Let u G p( m ) be quasi-definite and let {P n ) n >o be the associated sequence of 
left orthogonal matrix polynomials. Then, the following statements are equivalent: 

(i) u is semi-classical. 

(ii) There exist a polynomial a G P \ {0} with dga — p, s G N U {0}, s > p — 2 
and G C( m ' m ' (n > 0, -s < j < p), such that 

a(x)P^ +1 (x)= ^ n) Pn +j (x), 

j=-s 

where 6^ ^9 for some n> s (we use the convention P k = 9 for k < 0). 
Proof. 

(i) => (ii) If u is semi-classical then there exist a G P \ {0} and B G P^ m ' such 
that D(ual) = uB. Hence, (Pn+i) n>0 is quasi-orthogonal with respect to ual of 
order s —max{p — 2, q — 1}, where p =dga and q =dgB. This implies 

< x k P' 1 ,uaI >= 9, 0<k<n-s-l, 
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< x n ~ s P' n+l ,ual >^ 9, for some n>s. 
Since aP^ +1 G P^+p there exist 0^ verifying (ii). 

(ii) => (i) Let v = ual. Then 

p 

(x k P^ +1 ,v) = j n) (x k P n +j,u) = 9, 0<k<n-s-l, n>s, 
i=-s 

and, since < x n ~ s P n - Sl u > is non-singular, 

(x n - s P; i+1 ,v) =G^ s (x n - s P n - s ,u) for some n>s. 

Thus, (-Pn+i) n >o * s quasi-orthogonal 01 order s with respect to v. Obviously 
(P n ) n >o i s regularly quasi-orthogonal of order p with respect to v and, so, u is 
semi-classical, oo 

5.2 Theorem. (Diffcro-diffcrential equation) 

Let u G p( m ) be quasi-definite and let (Pn) n > be the corresponding sequence of 
left orthogonal matrix polynomials. Then, the following statements are equivalent: 

(i) u is semi-classical. 

(ii) There exist two polynomials a, (3 G P with p = dga > 0, q — dg/3 and 
matrices A ( k n) G C( m ' m ' (n > 0, -s < k < s) , such that 

a(x)P^(x)+(3(x)P^x) = J2 4 l) Pn + ^l 

k— — s 

where s >max{p — 2, q — 1} (we use the convention Pk = 9 for k < 0). 
Proof. 

(i) (ii) Let u be semiclassical. By Theorem 5.1, there exist a polynomial 
a G P, with dga = p\ > 0, a non-negative integer s\ > p\ — 2, and matrices 

£ C (m,m) ( n > Qj < j < ^ guch that 

Pi 

a(x)P' n (x)= ]T Q^Pn-i+jix), (5.1) 

j= — 8l 

Taking derivatives in (5.1), we obtain that 

a>P:: + aa>P> n = £ Q^aP^_ 1+J , 

j=Sl 

If we use (5.1) in the right hand side, it follows that 

pi pi 
a-P>; + aa>P> n = ]T £ P n -2 +j+k , 

j=—si fe=— Si 
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Let us denote a(x) = a 2 {x), (3{x) = a{x)a'{x), kf^ k = Q { .^ 1] Q { ^ 2+] \ and 
p = 2pi =dga, q = 2pi - 1 =dg/3, s =max{2pi - 2, 2si + 2}. Then, 



with s >max{p — 2, q — 1}. 

(m) =>■ (i) Let (7r„) >0 be the dual basis of (P n ) n >n , that is, 7r„ = uP^E^ 1 , 
where E n :=< P n ,P n >«• Let (Qn)„> be the basis of P( m ) given by Q n (x) — 

— ~^jPn+i> an( i l et {Pn)n>o be the corresponding dual basis. Thus, 

D Pn = -{n+l)n n+1 , n>0. (5.2) 
We consider for every n > the linear functional — D(n n a) + n n (3. There exists 

\ k Jk>0 

oo 

~D(ir n a) + ir n [3 = ^2 Pfc A i" } > 

k=0 

where from 



LXJ 

< Q J; -D(ir n a) + ir n f3 >= ^ < Q^p* > A^ = aJ 



k=Q 

holds for j > 0. Hence, our hypothesis implies that 



A ^ = ^TT E W,- (5.3) 

fc— — s 

For n = 0, (5.3) leads to A^ 0) = if j > s, and 

s-l 



k=0 

or 



s-l 



fe=0 



s-l 



-D(ua)P*E^ + u(3P*E^ = ^> fc A£ 0) . 
If we denote A^ 0) = A^E^P,)*)- 1 , then 

-£(ua)+ U /? = ^/> fe Af (5.4) 

fc=0 

In a similar way, for n = 1 (5.3) gives 

s 

-D(ua)P*E^ + uf3P*E^ = ]T p k X^ 



k=0 
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because = for j > s + 1. Thus, 



-D(uaP*) + uf3P* = Pk^k ] (5-5) 



with A^ = a£ 1} £7i. 

Keeping in mind that Pi (a;) = M\x + M2, with M\ non-singular, (5.5) remains 

s 

-D(ua)P* - uaM* + uf3P* = ^ p k \^\ 

k=0 

and, using (5.4), 

-uaM* = ± Pk X^-(Y /Pk X^)p*. 

fe=0 \fe=0 / 

Taking derivatives and applying (5.2) we obtain 
-D(ua)M* = -j2(k+l)n k+1 \W+[J2(k+ IK+iA^ j A*-(£pfeAf J = 

fc=0 \fe=0 / \fc=0 / ! 

= -E( fc + 1 )^Vi^ + \^ 1) +fE( fc + 1 )^Vi^+\^ 0) ) a*-(i>a£ 0)N ) mJ 

fe=0 \fe=0 / \fc=0 / 1 

As a consequence, the polynomial ^ G P^" 1 ^ given by 



*(a;) := 



_fc=0 \fe=o 
whose degree is not greater than s + 1, verifies that 

s-1 



D(ua) = 11$ + PfcA 
From this and (5.4) we get 



(0) 

k 

k=0 



D(ua) = uB 

* + 01 

with B = , that is, u is semi-classical. 00 



Remark 1. The distributional equation D(ua) = u/3 for a scalar functionals 
u implies that 

< x k {ap" n {x) +f3p' n (x)),u > = 
=< (x k p' n (x))' ,ua > -k < x k ~ 1 p' n (x),ua > + < x k p' n (x),uj3 >= 
= — k < x k ~ 1 p' n (x),ua >, 
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wich vanishes for k = 0, . . . , n — s — 1. So, 

s 

<XPni x ) + PPn( x ) = X nkPn+k(x) 

k— — s 

In the matricial case, the equality D(ual) — uB does not imply < x k P^(x), 
uB(x) >=< x k B(x)Pn(x),u > for the non-conmutativity, and more gencraly, 
there not exists a polynomial B(x;n) such that < x k P^(x),uB(x) > = 
< x k B(x;n)P' l (x),u >, neither. 

In this situation it is not possible to obtain a differo-differential equation with 
the polynomial B explicitly. 



5.2 Remark 2. Obviuosly, the differo-differential equation given by the The- 
orem 5.2 is not unique because we can modify it by adding any structure relation. 



§ 6 - Some examples. 

We will illustrate the preceding results with some examples. 

1. In the first one, we consider the scalar Laguerre weight w — xe^ x that verifies 
the Pearson type equation 

(wx)' — w{2 — x). 
Now, we define the matrix weight function 

that also verifies a Pearson type equation. In fact 

D(uxl) = uB(x) 



where B(x) 



2-x -2 
A~x f 
Moreover, there exists T G GL(R 2 ) such that 



u = TuT* 



xe~ x \ (1 

x 3 e- x r 



i.e., u is congruent with a diagonal weight with classical scalar weights in its diagonal 
and then, the functional u is positive definite. 

Let be the monic Laguerre polynomials asociated to the weight x a e~ x , a > 
-1. Then, 

'l ( n\x) 



Ln(x) - [ li?\x). 

constitutes a sequence of matrix O P related to the diagonal weight u. Obviously 
u satisfies a Pearson equation, 



D(uxl) = uB{x), B{x) = (T ) B{x)T 



2-x 
A-x 
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and a sequence of matrix P associated to u is given by 

P n {x) = L n {x)T, n>0 
It is easy to obtain a estructure relation 

xP^ +1 {x) = (n + l)P„+i + (n + 1) (" J 2 n ° 4 ) P„(x), 

as well as, for all n > 2, a differo-diffcrcntial equation, 

^'(x) + = n(n + l)P„(x) + 



„ + 3^»-'W'"r v ^ o (n + 2)(n + 3) 
Moreover, the polynomials satisfy the differential equation, 
xlP"(x) + (x-a-l)lP'(x)=rdP(x), 
and so, we can obtain for the polynomials P n the following differential equation: 

xK{x)+{^~ 2 x Q _^P' n {x)=nP n {x). 

For the monic polynomial P n = T~ 1 P n we have 

xPZ(x) +T" 1 ( x " 2 ^^TP^). „P„(x), 

and then 

a;P,;'(x)+ J B*(x)P; i (x)=nP n ( a; ) 



The polynomials P n can be obtained by means of a Rodrigues type formula 
Pn(x) = ^(T'T)" 1 u^ d -^ (3*T) . 



2. Now, we consider the the matrix weight u := e _:E R(x) where R(x) = 
j ^ ^ . So, is not positive definite. 

Then, u verifies the distributional equation 

v! = e- x2 R'(x) - 2xe- x2 R(x) = u(-2xl + R- 1 (x)R'(x)) = uB(x), (6.1) 
/ — 2x 1 \ 

where B(x) = ( ] . Moreover, u has not a congruent diagonal. 



Now, let 



P n (x) = (-2)(-")e- 2 5- 1 (x) d " (e 2l 5(x)) , n>0, (6.2) 
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1 + X 1 
x 1 
By derivation in (6.2) 



where S(x) 



P n+1 (x) = - l -[-2xI + S- 1 (x)S\x)]P n {x) l -P' n {x). 

Then, 

P n+1 (x) = ~[B*(x)P n {x) + P n (x)}, (6.3) 

and the leading coefficient of P n is / because Pq(x) = I. 

Applying the Leibnitz rule for the derivatives in (6.2), we have 

P n+2 ( X) = , 2: - V, :„^' : ^ ^ = 

= (-2)-^xe* 2 S-\x) dn+1 ^^ 

+ {-2)- n - 2 e x2 S- 1 {x)S'{x)S- 1 {x) i - 



_ lr d n+1 (e~ x S{x)) 



dx n + x 

So, we obtain the recurrence relation, 

P n +2(x) = -i[S*(x)P„ +1 (x) + (n + l)P n (x)]. (6.4) 

As a consequence, (P n ) is a sequence of matrix orthogonal polynomials with 
respect to certain non-singular functional (non-positive definite because —5(71+1)/ 
is a non- positive definite matrix). 

Moreover, from (6.3) and (6.4) we can obtain a structure relation 

P n {x)=nP n _ 1 {x), (6.5) 

and by derivation, the differo-differential equation 

PZ(x)=n(n-l)P n - 2 (x). (6.6) 

So, {P n ) is a semiclassical matrix orthogonal polynomial sequence. 
Now, from (6.3) and (6.5) it follows the differential equation, 

P n \x) + B\x)P^{x) = -2nP n (x). (6.7) 

Indeed, (P n ) is an sequence of matrix orthogonal polynomials with respect to 
the matricial weight u. In fact, 

/OO 
P' n {x)e- x2 R{x)P^{x)dx = 
-OO 

/OO 
K(x) + {-2xI + R'(x)R- 1 {x))P n (x)]e- x2 R{x)P t m {x)dx = 
-OO 
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/oo 
P n (x)e- x2 R{x)P^(x)di 
-OO 



the last equality from equation (6.7). In the same way, we obtain for the above 
integral that 

/OO pOC 
P' n {x)e- x2 R{x)P'^{x)dx = -2m / P n {x)e~ x ' '' R{x)P t m {x)dx. 
-oo J —oo 

So, when n ^ m we have the orthogonality for (P n ) with respect to u. 
As a consequence we have a quasi definite functional such that: 

(i) It satisfies the distributional equation D(ual) — uB. 

(ii) The corresponding sequence of matrix orthogonal polynomials satisfies a 
Rodrigues type formula. 

(iii) It verifies the structure relation (6.5), the differo-diffcrential equation (6.6) 
and a differential equation (6.7). Keeping in mind Remark 1 to Theorem 5.2, notice 
that we have (6.7) because uB = (uBf = B f u and B l P n = P n B l . 



3. In the last example, we consider the Jacobi scalar weight w = 1 — x 2 that 
satisfies the distributional equation 

(to(l - x 2 ))' = -w ■ Ax, x e[-l, 1]. 



We define the matricial weight function, 



u :- 



1 

x 



XG [-1,1], 



that verifies the Pearson type equation 

D(u(l-x 2 )l) =uB{x), 

where 

r>/ x ( -9x 2 + x 2 \ 

B ^ = { 1 -11* J- 

Notice that u is positive definite and then, there exists a sequence of matrix 
orthogonal polynomials (P n ) related to u. Moreover, u not is congruent with a dia- 
gonal weight. So, the sequence (P n ) satisfies a structure relation with s = 1, p = 2 
and a differo-differential equation of second order in the following way: 



2 

(1 - X 2 ) 2 P;:(X) - X(l - 2x)P' n { X ) = J2 ^nkPn+k- 

fe=-4 
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